Thermodynamic description of a dynamical glassy transition 
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For the dynamical glassy transition in the p-spin mean field 
spin glass model a thermodynamic description is given. The 
often considered marginal states are not the relevant ones for 
this purpose. This leads to consider a cooling experiment on 
exponential timescales, where lower states are accessed. The 
very slow configurational modes are at quasi-equilibrium at an 
effective temperature. A system independent law is derived 
that expresses their contribution to the specific heat. t/t w - 
scaling in the aging regime of two-time quantities is explained. 

25.10.Nr,25.40.Cx,25.50.Lk,64.70.Pf 



The structural glass transition occurs only in an ideal- 
ized adiabatic cooling procedure at the Kauzmann tem- 
perature Tk- In a realistic experiment with finite cool- 
ing rate a gradual freezing transition takes place in a 
small interval centered around a higher temperature Tf, 
that depends on the cooling rate. Although the freez- 
ing is not a sharp thermodynamic transition, there can 
be some 15 decades in time involved, from picoseconds 
to many hours. By extrapolation from the high and low 
temperature sides, one may define jumps in quantities 
such as the specific heat and the compressibility. It has 
been pointed out by Jackie [Q and Palmer Q that the 
freezing transition can be described as a smeared ther- 
modynamic transition, on a thermodynamic basis, with 
ensemble averages replacing time averages. The observa- 
tion time sets the time scale that separates "fast" pro- 
cesses (time scale less that the observation time) from 
"slow" ones. The latter are essentially frozen. 

Upon cooling, a liquid freezes dynamically in a glassy 
state with extensively smaller entropy. The free energy of 
the glassy state then is much larger, and it is not obvious 
why the system can get captured in such a state. The 
point is that the condensed system has lost the entropy of 
selecting one out of the many equivalent states; this part 
of the entropy is called the complexity, configurational 
entropy, or information entropy X. jlj] Q] This can be 
understood as follows. For long times the system is stuck 
in states with long (but finite) lifetime. These states are 
called "components" by Palmer Q and 'states' or 'TAP- 
states' in spin glass theory. When the Gibbs free energy 
F a of the relevant state a has a large degeneracy N a = 
exp(X a ), the partition sum yields Z = ^ a exp(— f3F a ) w 
N a exp(— ftFa), so F — F a — TX a is the total free energy 
of the system. The entropy loss arises when the system 
chooses the state to condense into, since from then on 



only that single state is observed. H As the total entropy 
S = S a + T a is continuous, so is the total free energy. 

The aim of the present work is to analyze this thermo- 
dynamical picture of the dynamical freezing transition in 
a well understood mean field model. The difficulty is to 
extract the information that is not due to the mean field 
limit. The gain is that we find strong constraints satisfied 
by the dynamics, without having solved the it. We first 
discuss the present status of marginal replica theory, its 
relation with dynamics and its fundamental flaws. Then 
we shall propose a solution to these paradoxes. 

We consider the mean field p-spin interaction spin glass 
model of N coupled spins in a field with Hamiltonian 
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The independent Gaussian random couplings J%^%2---i v 
have average zero and variance J 2 p\/2N P ~ 1 . The spins 
are subject to the spherical condition ^T Sf = N. 

This model has a close analogy with models for the 
structural glass transition. [Q On a static level there 
occurs a transition to a state with one step replica sym- 
metry breaking (1RSB) at a temperature Tk, comparable 
to the ideal glass transition in an adiabatic cooling exper- 
iment. The 1RSB replica calculation involves parameters 
three parameters. The overlap of spin configurations in 
two states can be equal the Edwards- Anderson parame- 
ter or self-overlap q±, or have the smaller value qo (= 
for H = 0); these values occur with probability 1 — x, 
and x, respectively. The free energy reads 



(2) 



where £ = 1 — x and Q = 1 — £,q\ — xqo ■ Unless stated dif- 
ferently, we shall take H = 0. qo and qi are determined 
by optimizing F. For x the situation is not unique, but 
depends on the time-scale considered. Setting dF/dx = 
yields the static phase transition at T = Tk- || When 
considering the Langevin dynamics of this model, one 
may derive dynamical equations for correlation and re- 
sponse functions by taking first first N — > oo. |7]] Solving 
for large t leads to a sharp phase transition at larger 
temperature T A = J{p(p - 2)P~ 2 /2(p - l)?" 1 } 1 / 2 , and a 
different form x(T). This dynamical value for x can be 
simply rederived from a replica calculation in which the 
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"replicon" or "ergodon" M fluctuation mode is taken 
to be massless. This leads to the marginality criterion 
p(p- l)(3 2 J 2 q\~ 2 (l-qi) 2 /2 = 1. At the dynamical tran- 
sition there is a sharp jump in the specific heat. |tJ 

To discuss the situation, we must first consider the 
TAP states. A state is labeled by a and has local mag- 
netizations mf = (Si) a . Its free energy F a (T) is a ther- 
modynamic potential that determines the internal energy 
and the entropy by its derivatives. In the present model 
F a = FxApi'm-i) is known explicitly. It is a minimum of 
the "TAP" free energy functional Q Q 



F TA p (mi) = - ^2 Jh- 



mi, ■ ■ - mi 



NT NBJ 2 
- — log(l - ? ) - -^(1 + (p l)<f - pq p ~ ) (3) 

where q — (1/^0 X)i m ? is the self-overlap. Below we 
shall argue that the commonly used Gibbs weight, p a — 
cxp(— (3F a (T))/Z is the relevant one. Given the PaS one 
can define the component averages such as F = J2 a PaF a , 
U = J2aPa U a, C = J2 a Pa C a, and even the complexity 
jjj Q X = — ^2 a p a l n Pa- For observables the direct eval- 
uation from the ordinary partition sum should coincide 
with the outcome of the TAP-analysis: U = U, M = M . 
They need not be derivatives of F . 

For all T < Ta the mode-coupling equations a massless 
ergodon, responsible for aging phenomena. We recently 
assumed that this is a very general phenomenon, will 
automatically get trapped in a state with diverging time 
scale, whenever present. The marginal replica free energy 
has the form 



F = F 



Tic 

x 



where 



V [ Ilog(p-l) + --l 
p 



(4) 



(marginality) (5) 



is the complexity of the marginal states. Below Ta the 
free energy lies below the one of the paramagnet and has 
a larger slope. This would naively imply a latent heat. 

There is, however, another prediction for the free en- 
ergy. |7]] It involves the internal energy and an entropy 
obtained by integrating (l/T)dU/dT from a temperature 
in the glassy phase up to some large temperature. The 
resulting "experimental" glassy free energy (?]] 
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exp 
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exceeds the paramagnetic free energy quadratically and 
is by construction a thermodynamic potential. It was 
reproduced by analysis of the TAP states. (l2) 

The difference between (|]) and (||) led us to ques- 
tion fundamentally the validity of replica calculations for 
dynamical 1RSB transitions. Our aim was to find the 



meaning of the logarithm of the dynamical replica free 
energy (eq. (||) with dF/dqo — dF/dqi — but with 
dF/dx^O). 

By doing the full thermodynamic analysis of the TAP- 
partition sum ^ a exp(— f3F a ) at H = 0, we found that 
that the replica free energy ([!]) is reproduced. |TJ] As 
the glassy free energy lies below the continuation of the 
paramagnetic one, we considered this as proof that the 
complexity is the driving force for the dynamical phase 
transition. [Q In the doing the analysis we realized that 
the calculation of the dynamical complexity cannot be 
separated from the calculation of the free energy, and 
that replica symmetry breaking is essential. The prob- 
lem boils down to a replicated TAP free energy that has 
6 replicated order parameters. We have now extended 
this analysis to H > 0. For 1RSB with a common break- 
point x, each replica order parameter now brings 3 pa- 
rameters. We thus obtain an optimization problem in 
18 variables, that was solved partly using Maple. We 
have verified that the total free energy, the internal en- 
ergy and the magnetization, calculated within the TAP- 
approach, coincide with their replica values. For the 
magnetization this is particularly satisfying, as M a of 
a given marginal TAP-state is temperature independent. 
(It then holds that 8M a /dT ^dS a /dH = 0. Neverthe- 
less, the component average M — ^2 a p a M a is tempera- 
ture and field dependent, and equal to the replica value 
M = pH{l - - xq \y 

For interpreting eq. (^) one might be tempted to con- 
sider I c /x as the full complexity. Since x — > for T — > 0, 
this is hard to explain on a physical basis, however. The 
1 jx dependence in (^) does not disappear after the quan- 
tization of the spherical model, recently proposed by 
us. p8pl Analysis of the equations for quantized spher- 



ical spins, or for Ising spins, learns us that, though Ta 
shifts, the term TI c /x(T) survives for marginal states. 

The other interpretation of eq. (^) is that T e =T/x is 
an effective temperature at which the slow processes lead- 
ing to I c are in quasi-equilibrium. This interpretation is 
promising, since the same effective temperature shows up 
in the fluctuation-dissipation relation in the aging regime 
of the mode-coupling equations. M 

There remain some paradoxes connected to the 
marginal states. For large p one has I c ~ (N/2)\ogp, 
which (for quantized spherical spins or for Ising spins) 
already exceeds the total entropy available. This shows 
that the dynamical transition at Ta has no thermody- 
namic counterpart in short range systems, at least for 
large enough p. In fact, Ta may be identified with the 
critical temperature of mode-coupling theory, which lies 
well above the freezing temperature TV. Even more cum- 
bersome is that for marginal states with H ^ we have 
shown violation of the inequality C > C Q. This says 
that marginal states are intermediate-time states, from 
which the system must escape at longer times. 

The implication of these arguments is rather dramatic: 
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for comparing with realistic short range systems we must 
abandon the assumption that the marginal states are the 
physically relevant states. Let us see how this could hap- 
pen. Marginality arises automatically in the dynamical 
equations after taking first the limit N ~ * oo and then 
t — ► oo. This order of limits, however, prevents all ac- 
tivated processes, which in the mean field model would 
need a time ~ exp(TV). As a result, all dynamics is con- 
fined to the highest TAP-states, which are marginal. The 
lower states are never reached. In a realistic short range 
glassy system, however, there is no sharp distinction be- 
tween slow and activated processes, and the latter can 
certainly not be omitted. So in a realistic experiment we 
do expect to reach lower states. 

In order to compare to realistic systems and to avoid 
thermodynamic paradoxes, we propose another look at 
the mean field system. We consider the system at fixed 
large N under such conditions that a range of TAP-states 
below the highest (marginal) ones are accessed. We thus 
focus our attention on the free energy of the state, and 
not so much on the timescale needed to reach it. At H = 
the free energy of the TAP states can be characterized 
by a parameter 77 (rj s t < 77 < 1), that enters the condition 



P 2 P {p-l)qV- 2 (l-q) 2 /2 = n. 



(7) 



For r\ = 1 this would be the marginality condition, while 
the static equations can also be put in this form j5j with 
V — Vst < 1 independent of T. One now finds the break- 
point x = (1 — q)(p — 1 — T))/qrj, so choosing 77 between 
rj st and 1 can alternatively be seen as a way of fixing 
the mysterious parameter x, for which no obvious crite- 
rion was present. In the present approach it is directly 
related to the time scale t v = exp(Nr v ) at which the 
?7-states are reached. Unfortunately, the precise relation 
between the logarithmic time variable t v and the lowest 
reachable TAP free energy at that scale, F n , is unknown. 

We assume that at a given exponential timescale t v the 
dominant states are determined by a saddle point. Such 
a behavior was seen in a solvable glassy transition in a 
directed polymer model, recently proposed by us. |2l]] 
This allows to restrict the discussion to TAP-states with 
a common free energy. 

We have verified that the equivalence between the 
replica and the TAP analysis also holds for 77 < 1. The 
77-states have complexity 

^)-^£zl- feziz^±a } ( 8 ) 

2 77 pr\ 

As r\ decreases from 1 to n st , 1 c {rf) goes down from the 
marginal value I c (l) to the static value T c (n st ) = 0. A 
sensible cooling experiment must have a temperature de- 
pendent value of 77, such that 77 — > 77 st fast enough for 
T — > 0. In hindsight this is just what is needed to de- 
scribe a realistic cooling experiment in the mean field p- 
spin model. Let us assume to have some logarithmically 



slow cooling traject T(t) where also 77 = 77(f) depends on 
time in a still unknown, but determined fashion. We can 
then eliminate f and construct the function T(rj) (and its 
inverse r](T)) that characterizes our cooling traject. A 
dynamical freezing transition will occur when the lowest 
state reached at time f freezes at the temperature T(t). 
This occurs at Tf = T(i~if) where 77/ is the solution of 



T(r,) = j 



p 7 i{p ~ 1 - v) p 

2(p - I)?" 1 



2 1 V2 



(9) 



As T(t) and 77(f) will depend on the cooling procedure, 
we expect T{rf) to do the same, so that Tt will not be 
universal but depend on the specific traject. 

Also the assumption of Gibbs weights is now justi- 
fyable. TAP states with free energy larger than or equal 
to the ones fixed by 77 are now effectively in thermody- 
namic equilibrium, and may be described by the Gibbs 
weight. Lower ones play no role anyhow. 

The paradoxes related to marginal states, signaled 
here, occur much more widely. At present it is a whole 
field of research to consider dynamics of mean-field mod- 
els by first taking the mean field limit and then consid- 
ering large times. Though the approach has relevance 
for short and intermediate time dynamics, its long time 
regime is a result of "squeezing the system into marginal 
states" , which has no bearing on the long time relax- 
ation of short range models. This is already expressed 
by the unique sharp dynamical transition temperature 
Ta, found in the dynamical approach. It disappears on 
exponential time scales, and is replaced by a cooling-rate 
dependent freezing temperature. 

The basic object in our approach is the complexity. An 
important question is whether it can be measured in the 
glassy phase. If so, it should be related to the specific 
heat or the temporal energy fluctuations. When moni- 
toring the internal energy as function of time, as is easily 
done in a numerical experiment, one obtains essentially a 
noisy telegraph signal. Each plateau describes trapping 
in a TAP-state for some definite time. The variance of 
the noise in the internal energy on this plateau is equal to 
T 2 C=T 2 J2p a C a =J2 a ((6U a ) 2 )=T 2 j:p a dU a /dT. From 
time to time the system moves to another TAP-state, 
causing additional noise. The variance of the total noise 
equals T 2 C, and it should exceed T 2 C. @ In ref. @ 
it was pointed out that there can be an extensive differ- 
ence between the specific heat C = dU/dT = dU /dT = 
J2 a d(p a U a )/dT and the component average energy fluc- 
tuations C—TdS/dT We can now consider their differ- 
ence at H = in a cooling experiment of the type intro- 
duced above. We find the excess specific heat 



AC = C -C = N- 



3(1 - v) d V 



(10) 



2pr)(l-q) dT 
Using this we can easily verify the fundamental relation 
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^1 C (T) = ±x(T)AC(T) 



(11) 



It expresses the complexity in terms of measurable quan- 
tities, namely the excess specific heat and the ergodicity- 
breaking parameter x. Equation ( pi] ) holds equally well 
for H ^ but fixed. In the generalization of ( |To| ) one now 
encounters parameters ry(T; H) and its derivative drj/dT. 
The proof of eq. ( plf ) for H ^ is lengthy but could be 
verified using Maple. [jl6| 

The complexity can also be measured along the tran- 
sition line in the (T, iJ)-plane. This issue is related 
with the Ehrenfest relations, and has been discussed else- 
where. JT|] 

With T e = T/x eq. (O) can also be written as 



(12) 



dU dS dl c 
dT dT e dT 



Our interpretation of the replica results leads to an 
effective temperature T e (t) — T(t)/x(t). The slowest ac- 
tive modes are at quasi-equilibrium at this effective tem- 
perature. This explains why T e also shows up in the 
fluctuation-dissipation relation. g|, jl5| As they set the 
slowest time-scale, they must also dominate the dynam- 
ical free energy. This is why the change of structural 
modes dX c /dT in eq. ( p^ ) has prefactor T e . 

In numerics on the fluctuation-dissipation relation in 
spin glasses and even binary soft spheres ji~7| it has been 
observed that the factor x is linear in T . Let us give 
a simple explanation for that behavior. As the effective 
temperature must exceed the Kauzmann temperature, 
we can estimate T e (T) const, implying indeed x = 
T/T e ~ T. 

If we quench the system from high temperatures and 
let it evolve freely during a waiting time t w , T e will be 
set by equating the equilibrium relaxation time r eq at 



the instantaneous T P to t v 



' eq 



(T e ) — t w . Naively, one 



expects two-time quantities in dynamics to be a function 
of (t — t w )/T eq (T e ) = t/t w — l, explaining immediately the 
often observed t/t w scaling in the aging regime. There 
could be logarithmic corrections to this behavior. 
Replacing x — > T/T e we can write eqs (0) and (0) as 



F(T,T e 



U -TS- T e X c 



(13) 



It is a dynamical free energy determining S = —dF/dT 
and l c = -dF/dT e , while U = F + TS + T e l c . The 
system- independent laws (12) and ( |l3| ) are the corner- 
stone for our thermodynamic description of the dynami- 
cal glassy transition, and expected to be valid in general. 

We have verified eq. ([l2|) for a cooling procedure 
in the Ising chain with Glauber dynamics. At T = 
it is a coarsening problem of alternating up and down 
clusters of average length £(t) = V '4mt and energy 
U(t) = NJ(—1 + 2/£). counted by standard manners. 
Its logarithm reads T c — N(l + ln£)/£. From the inter- 
nal energy one may introduce an effective temperature: 



l/£ = e - 2l3 ' J ^ T e = 2J/ln£(t). The same value would 
follow when defining T e from the complexity T c . Both re- 
sults can thus be combined in terms of a dynamical free 
energy F = — iVT e m(2cosh/3 e J). This is a special case 
of eq. @ with 5 = 0. 

As above, we then consider a cooling experiment where 
T(t) = 2x(t)J/ln£(t) with a smooth decreasing function 
x(t). In the initial regime x(t) > 1 the system will achieve 
equilibrium at the instantaneous temperature T(t). For 
x(t) < 1 this will not happen. The system falls out of 
equilibrium and behaves as at T = 0: it is at quasi- 
equilibrium at T e (t) — T(t)/x(t). The freezing transition 
occurs around Tj = T(tf) where x(tf) = 1. In the frozen 
phase eqs. (||Jl§) are valid with S = 0. At T = (and 
for all T < Tf) the on-site correlation function is found 
to be a scaling function of t/t w , 

Summarizing, we have shown that a thermodynamic 
description of a dynamical freezing transition can be 
given. IH We have been led to discard the whole issue 
of marginal states. For comparing with realistic short 
range systems, we consider the mean field system at ex- 
ponential time scales, where lower states are accessed. 
They contribute to the partition sum of Gibbs weights 
over dynamically relevant states. This approach natu- 
rally leads to slow cooling procedures where a dynamical 
freezing transition occurs at a tunable temperature. This 
dynamical transition is described by a free energy that 
depends on the real temperature and on an effective tem- 
perature. 
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